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1. Introduction. 

Several groups are working on the development of a new generation of frequency and time 
standards based on atomic optical transitions [TJ [21 Neutral ytterbium atoms loaded 
into a laser lattice are very promising candidates for constructing such high-performance 
atomic clockjU [5]. The properties of Yb2 are necessary to assess the feasibility of using 
laser cooled and trapped Yb atomic species for ultraprecise optical clocks or quantum 
information-processing devices [6]. 

Unfortunately, reliable experimental data on the dissociation energy, equilibrium 
distance, and spectroscopic constants of the Yb2 molecule are unknown; e.g., the 
uncertainty of the experimental dissociation energy estimate from [7], 0.17 eV, is 
comparable to the value itself. A series of papers was devoted to their calculation. 
Dolg and co-workers reported the values 400 [8], 470 [9], and 740 cm -1 [10] for the 
dissociation energy. The coupled electron pair approximation, density functional theory 
(DFT), configuration interaction method with single and double excitations (CISD), and 
coupled cluster method with single, double (and non-iterative triple) cluster amplitudes, 
CCSD(T), were used to account for the correlation effects. Ytterbium core shells were 
replaced by scalar (spin-averaged) energy-consistent pseudopotentials (PPs) generated 
for 2, 10, and 42 explicitly treated electrons. Their latest result [10] of 740 cm -1 should 
be considered the most reliable, because the basis set used was larger than that in the 
previous calculations and 42 electrons were treated explicitly for each Yb atom (42e-PP). 

In [UJ, the scalar DFT approach and a 24-electron relativistic effective core 
potential (RECP) model for ytterbium were used. The dissociation energy estimates 
ranging from 500 to 1400 cm -1 were obtained with different exchange-correlation 
functionals. In [12], the Yb dimer was studied within the averaged quadratic coupled 
cluster, CCSD(T), and DFT approximations. The Yb atom was described by a 42- 
electron energy-consistent PP. It was emphasized that the "incomplete convergency, 
most clearly seen for Yb2 results, indicated the need for more advanced ab initio 
schemes". In addition to the evident problem of the incompleteness of the one-electron 
basis set, it is not clear whether the truncation of the cluster expansion after the three- 
body terms provides a good approximation for the Yb2 ground state, which could be 
considered as a perturbed four-electron system. Furthermore, almost all calculations 
mentioned above were done within the scalar relativistic approximation. In spite of the 
closed-shell-like nature of the system under study, the contribution of spin-dependent 
interactions to the bond energy can still be significant (cf. [13]). The only attempt [S] 
to estimate the role of spin-orbit interactions in Yb2 was made within a somewhat 
simplistic four-electron CI scheme using a very restricted basis set that can hardly be 
used to reliably reproduce the van-der-Waals behaviour of the potential curve. 

In the present paper, we report our results of improving the accuracy of the 
calculated dissociation energy, equilibrium distance and spectroscopic constants for 
Yb2 using extremely flexible generalized correlation basis sets, contributions from high- 
order cluster amplitudes and spin-dependent relativistic effects. Such improvements 
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were successfully applied to accurate calculation of Hg 2 when giving a few times better 
agreement with the experimental data than the other earlier performed studies [14]. 

2. Calculations and discussion. 

Scalar relativistic calculations were performed within the generalized relativistic effective 
core potential (GRECP) model[I5|[T6l[T71[Tg] using the CCSD(T) method (implemented 
in the molcas program package [19]) for correlation treatment. The high accuracy and 
reliability of this approach has been demonstrated in similar calculations (see, e.g., 
[T3]). The 4/ 5 /2, 4/7/2, and 6si/ 2 spinors of the Yb atom have the one-electron energies 
of -0.54, -0.49, and -0.20 a.u., respectively, and are usually considered as valence ones. 
The average radii of the 4s 1/2, 4pi/ 2 , 4p 3 / 2 , 4<i 3 / 2 , and 4rf 5 / 2 spinors are very close to 
those of the 4/5/2 and 4/7/2 spinors. Therefore, excitations of electrons from the latter 
spinors will lead to strong relaxation of the former spinors which, therefore, should be 
considered as the outercore ones for the GRECP generation procedure and "valence- 
type" property calculations[15j. Thus, we use the GRECP with 42 explicitly treated 
electrons for each Yb atom. In series of preliminary calculations, we have estimated 
the contributions from correlations with different shells of Yb to the dissociation 
energy of Yb2 (some of them are presented in table [3]). The main contribution is 
provided by the 6s shell whereas the contributions from the 4/, 5s, and 5p shells are 
relavitely small. It is clear that the corresponding contributions from the innermore 
4d, 4p, 4s, etc. shells will be significantly smaller. Thus, we "freeze" the 4s, 4p, 4d 
or 4s, 4p, id, 4/, 5s, 5p shells in 48- and 4-electron calculations, respectively, to reduce 
computational expenses in the present calculations. Generalized correlation basis sets 
comprising (19, 17, 7, 17, 6, l)/[7, 8, 5, 4, 3, 1] functions, basis C (core), in the former and 
(38,22,24,14,7,1)/[11,10,9,7,5,1], basis L (large), (38, 22, 24, 14)/[5, 6, 4, 3], basis M 
(medium, with the g and h harmonics removed from the previous uncontracted basis 
set), (38, 22, 24, 14)/[5, 5, 3, 2], basis S (small), in the latter cases were generated by the 
procedure developed previously [20l |2T] . 

Calculations were carried out for internuclear distances from 6 to 14 a.u. All our 
results were rectified using the counterpoise corrections (CPC)|22], [23] calculated for 
the Yb 6s 2 state with one more Yb atom treated as the ghost one. The energies of the 
rovibrational levels were obtained by solving the rovibrational Schrodinger equation with 
the numerical interatomic potential by the second order finite-difference method|24j. 
The stage of calculation of the molecular constant [21] begins with fitting the numerical 
potential curve for the dimer by polynomials with the help of the quasi-Hermitian 
method. Appropriate derivates of the potential curve at the equilibrium point are 
calculated by recurrence relations. Then rovibrational Schrodinger equation is solved 
by the Dunham method to express the Dunham coefficients in terms of these derivates. 

The X E+ closed-shell ground state of the Yb2 molecule disscociates into two Yb 
atoms in the 4/ 14 6s 2 ( 1 5') ground state. The computed ground-state potential energy 
curves for the Yb2 molecule are shown in table [1] and figure [lj the energies of the 
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lowest rovibrational levels for the ground electronic state are collected in table [2J our 
estimates for the dissociation energy, equilibrium distance, and main spectroscopic 
constants are listed in table [3j We started from 4-electron scalar relativistic CCSD(T) 
(denoted as 4e-CCSD(T) below) calculations with rather large basis set L, which gave 
D e =706 cm -1 . The negligible CPC (0.3 cm -1 for dissociation energy) indicates a good 
quality of the basis set used. Subsequent calculations of the effects of the difference 
between the iterative and non- iterative triple cluster amplitudes (CCSDT-CCSD(T) or 
contribution from iteration of triples) as well as of quadruple cluster amplitudes (these 
two contributions are denoted further as the iTQ contribution), valence - outer core 
correlations (OC), and the spin-orbit interaction (SO) described below have shown that 
the corresponding contributions to the Yb-Yb interaction energy are within 15% (with 
respect to the final dissociation energy estimate of 786 cm -1 ), thus justifying the choice 
of the 4e-CCSD(T) scheme as a good initial approximation. Note that the 4-electron 
FCI or 48-electron CCSD(T) calculations with considerably smaller basis sets M or C 
have given much lower D e estimates (D e =536 or 353 cm -1 , correspondingly). Thus, the 
quality of the basis set is of crucial importance for accurate calculations of the ytterbium 
dimer. 

The contribution from the quadruple cluster amplitudes as well as the difference 
between the iterative and non-iterative triple amplitudes was estimated as the difference 
between the total energies obtained in the 4e-FCI and 4e-CCSD(T) calculations with 
basis set M for each of the above mentioned internuclear distances. This difference was 
then added to the total energy obtained in the 4e-CCSD(T) calculation with basis set L. 
The derived correction from the iTQ amplitudes to the dissociation energy, 117 cm -1 , 
is 15% with respect to our final value for D e . We expect that the contribution from the 
iTQ amplitudes calculated with basis set M will change very slightly (with respect to 
the final D e value, etc.) if this basis set is replaced by basis L. It should be noted that 
the 4e-CCSD energy difference for R e = 100. a.u. and R e = 9. a.u. (the latter is close to 
the equilibrium distance) changes by 229 and 215 cm -1 in going from basis set S to M 
and from basis set M to L, respectively. The corresponding contributions from the non- 
iterative triple cluster amplitudes are 94 and 81 cm -1 , whereas the former contribution 
from the iTQ amplitudes is only 18 cm -1 . Thus, extrapolation to the infinite basis set 
limit should only slightly increase the dissociation energy estimate. 

The contribution from the correlations with the 4/, 5s, and 5p outer-core electrons 
was estimated as the difference between the total energies found in the 48e-CCSD(T) 
and 4e-CCSD(T) calculations with basis set C for each of the above mentioned 
internuclear distances. The only difference between these two calculations is the 
number of correlated electrons, therefore, the lowerings in the total energies give the 
contribution of the OC correlations. These lowerings were then added to the 4e- 
CCSD(T) and 4e-CCSD(T)+iTQ total energies derived above. The 4e-CCSD(T)+OC 
and 4e-CCSD(T)+iTQ+OC dissociation energy, equilibrium distance and spectroscopic 
constants were calculated with the obtained potential curves. The dissociation energy 
was decreased by 56 cm -1 , whereas the corresponding CPC contribution was about 
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of 100 cm -1 . It should be noted that the level of the approximations made in 
calculations [TOl [12] most closely correspond to our 4e-CCSD(T)+OC approximation 
(see table [3]). The difference (of 82-98 cm" 1 for D e ) between the results of the above 
calculations can be assigned to the insufficient flexibility of the basis set functions from 
the "outer-core" region in [TOl [12] . Unfortunately, calculations for 4 correlated electrons, 
which are necessary to check this conjecture, were not reported in the cited works. 

The effect of spin-dependent (effective spin-orbit) interactions was taken from 
Ref. [25] as the difference between the ground-state potential curves obtained in two- 
component relativistic DFT calculations |26j with full RECPs and in scalar relativistic 
DFT calculations with spin-averaged RECPs. The details of the employed procedure 
can be found elsewhere |2 7]. We only note that effects of electronic correlations are 
taken into account within DFT, so the additivity of the correlations and spin-orbit 
effects is irrelevant to (is not exploited in) our present study. An uncontracted Gaussian 
basis set [25] (10sllp8<i9/4g) was used to expand auxiliary one-electron spinors in the 
Kohn-Sham scheme. Two generalized gradient approximations for exchange-correlation 
functionals were employed, a rather universal Perdew-Burke-Erzernhof (PBE) model [28J 
and the Perdew-Wang approximation (PW91|29J), which is often believed to be 
particularly well suited for the description of van-der-Waals bonds [30]. 

Scalar relativitic DFT dissociation energy is about 1.6 times higher than that of the 
corresponding ab initio 4e-CCSD(T)+iTQ+OC approximation. Nevertheless, we believe 
that the spin-orbit contribution is rather correctly extracted from the DFT calculations 
because the spin-orbit interaction is described by the the one-electron operator and the 
one-electron parts are the same for Kohn-Sham and Schrodinger (including Hartree- 
Fock) Hamiltonians. Despite the potential energy functions obtained in the relativistic 
DFT calculations with PBE and PW91 functionals are slightly different in shape 
(dissociation energy estimates in two-component calculations are 1238 and 1298 cm -1 , 
respectively), the corresponding spin-orbit corrections to bond energies as functions of 
the internuclear separation almost coincide. The addition of these corrections to the 
results of accurate scalar relativistic calculations has increased D e by 19 cm -1 . 

3. Conclusions. 

We predict the exact D e and w e to be slightly higher than 786 cm -1 and 24.1 cm -1 and 
the exact R e to be slightly lower than 4.582 A , because all contributions (taken 
into account in our calculations with a good accuracy) except for the OC correlations 
change these constants corresponingly. We expect that the reported estimates of the 
dissociation energy, equilibrium distance and spectroscopic constants of Yb2 obtained 
by the CCSD(T) technique with very extensive basis sets and the incorporation of 
corrections for higher-order cluster amplitudes and spin-orbit interactions are the most 
reliable up to date. Our analysis has revealed a non-negligible role of quadruple 
amplitudes as well as the significant contribution from iteration of triple amplitudes 
in the cluster expansion (which were not taken into account in [8j [9j UQl U2] ) and small 
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but non-negligible contributions from spin-dependent relativistic effects. 
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Table 1. Potential energy functions for the Yb 2 ground state calculated with the help 
of the GRECP and different correlation methods, fntcrnuclear distances R and total 
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Table 2. The energies of the lowest rovibrational levels for the 1 E+ ground state of 
the 171 Yb 2 molecule from the GRECP/4c-CCSD(T)+iTQ+OC+SO potential curve 
are in cm -1 . J and v are the rotational and vibrational quantum numbers. 
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Table 3. The dissociation energy, equilibrium distance, and spectroscopic constants 
of the 1 E+ ground state of the 171 Yb2 molecule. R e is in A, B e in 1CT 3 cm , a e in 
10 -5 cm -1 , Yq2 in 10~ 9 cm -1 , and other values in cm -1 . 
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Figure 1. Calculated potential energy functions for the Yb2 ground state. 
Curve (a) corresponding to the computational scheme 4e-CCSD(T)+OC provides 
an approximation for all-electron scalar CCSD(T), that obtained at the 4e- 
CCSD(T)+OC+iTQ level (b) should approach the scalar relativistic limit whereas 
curve (c) presents our best full relativistic results (4e-CCSD(T)+OC+iTQ+SO). 



